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Context. A photoionized gas in thermal equilibrium can display a thermal instability, with three or more solutions in the multi- 
branch region of the S-shape curve giving the temperature versus the radiation-to-gas-pressure ratio. Many studies have been 
devoted to this curve and to its dependence on different parameters, always in the optically thin case. 
rS^ • Aims. The purpose of our study is the thermal instability in optically thick, stratified media, in total pressure equilibrium. We 
are also interested in comparing photoionization models issued from the hot and cold stable solutions, with the currently used 
^ , models computed with an approximate, intermediate solution. 

5^ ' Methods. We have developped a new algorithm to select the hot/cold stable solution, and thereof to compute a fully consistent 
photoionization model. We have implemented it in the TITAN code and computed a set of models encompassing the range of 
03 conditions valid for the Warm Absorber in Active Galactic Nuclei. 

Results. We have demonstrated that the thermal instability problem is quite different in thin or thick media. Models computed 
. with the hot/cold stable solutions, or with an intermediate solution, differ all along the gas slab, with the spectral distribution 
■ changing as the radiation progresses inside the ionized gas. These effects depend on the thickess of the medium and on its 
5—j ' ionization. 

Conclusions. This has observational implications in the emitted/absorbed spectra, ionization states, and variability. However im- 
possible to know what solution the plasma will adopt when attaining the multi-solutions regime, we expect the emitted/absorbed 
spectrum to be intermediate between those resulting from pure cold and hot models; such a phase-mixed medium can be well 
reproduced by intermediate solution models. Large spectral fluctuations corresponding to the onset of a cold/hot solution 
could be observed in timescales of the order of the dynamical time. A strong turbulence implying supersonic velocities should 
permanently exist in the multi-branch region of thick, stratified, pressure equilibrium media. 

Key words. Instabilities - Radiation mechanisms: thermal - Radiative transfer - Methods: numerical - Galaxies: active - X-rays: 
general 



1. Introduction 

It is well known that a photoionized gas in thermal equilib- 
rium, i.e. a gas where the radiative heating is balanced by 
radiative cooling, can display a thermal instability (e.g., 
Krolik, Tarter & McKee 1981). The phenomenon mani- 
fests itself in the S-shape of the net cooling function^ or, 
which is equivalent, through the curve giving the temper- 
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1 In this study, the net cooling function L is defined as the 
difference between the cooling (A) and heating (r) functions 
in a given region, and is expressed in units of erg cm 3 s _1 . It is 



ature versus the radiation-to-gas-pressure ratio. Such an 
S-shape curve allows for the co-existence of gas at different 
temperatures and densities for the same pressure ratio. 

For a given value of the radiation-to-gas-pressure ra- 
tio, the gas can then be in three (or even more) states of 
thermal equilibrium, which depend on the Spectral Energy 
Distribution (SED) of the ionizing spectrum, on the abun- 
dances, and on other physical parameters affecting the 
heating and the cooling. One of these states is thermally 



also possible to use a net cooling function describing the energy 
balance per gram of material, per second. 
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unstable, as it does not satisfy the stability criterium to 
isobaric perturbations (Field 1965): 

BC 

(^)p gaa > 0, (1) 

where C is the net cooling function, T is the temperature, 
and P gas is the gas pressure. The two other states for which 
the derivative of C is positive, are stable. Such thermal 
instability is, for instance, at the origin of the two-phase 
model for the interstellar medium, where cold (T < 300 K) 
neutral atomic and molecular clouds are embedded in a 
warm (T ~ 10 4 K) intercloud medium (Field, Goldsmith 
& Habing 1969). 

1.1. Optically thin media 

Krolik et al. (1981) studied the thermal instability in the 
context of Active Galactic Nuclei (AGN), for optically 
thin, X-ray illuminated gas. In such media, the radiation 
pressure keeps a constant value, and the energy balance 
equation solution depends on the gas pressure only. When 
the gas pressure is small compared to the radiation pres- 
sure, there is a unique, stable, "hot" solution, where both 
the heating and the cooling are dominated by Compton 
processes (Compton heating, inverse Compton cooling). 
As the gas pressure increases, atomic processes (photoion- 
ization heating, line and continuum cooling) become im- 
portant, and multiple solutions arise. The harder the spec- 
trum, the more extended is the region encompassed by the 
multiple solutions. Then, above a given gas pressure, again 
a unique, but this time "cold" solution, arises. Krolik et 
al. showed that such multiple solutions may lead to the ex- 
istence of a two-phase medium consisting of a hot, dilute 
gas confining denser, cooler clumps, which they identified 
as the Broad Line Region in AGN. 

To better address the thermal instability issue, these 
authors introduced the ionization parameter 2 (dimen- 
sionless), defined as 



where F- lon is the incident ionizing flux between 1 and 
10 3 Rydberg, n is the total numerical density, and T the 
gas temperature. In a fully ionized gas, the radiation-to- 
gas-pressure ratio, P r ad/-Pgas, is about equal to 2/2.3. The 
so-called S-shape curve is thus, actually, the curve giving 
T versus 5. Krolik et al. showed that the multi-branch 
region of this curve (or, which is equivalent, the multi- 
solutions regime) corresponds to a range of 5 between 1 
and 10 for a typical AGN spectrum. 

Since then, many studies have been devoted to the S- 
shape curve and to its dependence on the incident ionizing 
spectrum and on the abundances, but always in the opti- 
cally thin case. The purpose of our study is different from 
previous works, as we are interested in optically thic^ me- 
dia such as the Warm Absorber (WA) in type 1 AGN, the 

2 Note that by "optically thick", we here mean media opti- 
cally thick to the photoionization continuum; this is the case 



X-ray line-emitting gas in type 2 AGN, or irradiated accre- 
tion discs in AGN and X-ray binaries. Our studies may be 
applied to media in any pressure equilibrium conditions, 
e.g., constant gas pressure, constant totajfl pressure, or 
hydrostatic pressure equilibrium. 

1.2. Optically thick media 

Let us assume a medium in total pressure equilibrium, 
consisting of a slab of gas with a given density at the illu- 
minated surface. Under such conditions, irradiation by X- 
rays induces a thermal instability beyond a given layer in 
the slab. The shape of the radiation spectrum and the ra- 
diation pressure itself, depend on the considered position 
in the gas slab. Near the surface, the radiation spectrum is 
close to the incident spectrum irradiating the slab (close, 
yet different, owing to the additional presence of radiation 
returned by the non-illuminated side of the slab) . 

If the radiation-to-gas-pressure ratio is sufficiently 
large, the equilibrium temperature is high, and these lay- 
ers are hot, displaying a unique, "hot" solution. Such equi- 
librium temperature stays almost constant in this region. 
In deeper layers, two effects take place. First, the radiation 
spectrum is harder, owing to soft X-ray absorption by the 
previous layers; it contains only hard X-rays and diffuse 
ultraviolet radiation created within the slab. As a conse- 
quence, the S-shape curve is more pronounced. Second, 
in the conditions of our models (described in Sec. [3]), the 
radiation flux decreases — and so does the radiation pres- 
sure — entraining the gas to enter the multi-solutions 
regime; the gas has now the choice between two stable 
solutions, a "hot" and a "cold" one. Finally, as the radia- 
tion penetrates still deeper in the slab, it suffers more and 
more modifications, and we get to a point where there is 
again a unique, "cold" solution. The corresponding equi- 
librium temperature remains almost constant again. This 
phenomenon can occur several times in the slab, accord- 
ing to the spectral distribution of the specific intensity 
through the medium. Each time the gas-phase changes, 
adjusting to a colder solution, the temperature profile dis- 
plays a sharp drop by one order of magnitude or so. In 
between phase changes, the temperature remains more or 
less constant. 

1.3. Electronic conductivity and alternative methods 

One of the questions one may ask at this stage is "which 
solution to choose when there are two or more stable 
ones?" The answer depends on the previous history of 



for gas irradiated by a continuum that gets absorbed at one 
or more ionization edges, being therefore altered when passing 
through the medium. In a weakly ionized medium, this can 
occur for a column density as small as 10 18 cm~ 2 . 

3 In this paper, the total pressure includes the gas and radi- 
ation pressure, only; however, it would be possible to include 
other contributions, like a turbulent and/or a magnetic pres- 
sure component to the total pressure. 
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the medium. A unique equilibrium state could be found 
if electron conduction was to be included (Begelman & 
McKee 1990; McKee & Begelman 1990; Rozahska 1999). 
The well-known numerical difficulty linked to the ther- 
mal instability could be overcomed using an integral for- 
mulation of the criterium giving the equilibrium solution, 
instead of a differential one (Rozahska & Czerny 2000). 
However, coupling conductivity effects with a complete 
radiative transfer in photoionization models is a compli- 
cated task, and has not yet been performed. Attempts to 
introduce elcctonic conduction are presently being under- 
taken by Chevallier et al. (in preparation). 

There are, however, other means to circumvent the 
problem of multiple stable solutions when dealing with 
photoionization codes. Two possibilities are discussed 
here: the first one was used to compute the vertical struc- 
ture of irradiated accretion discs in hydrostatic equilib- 
rium (e.g., Ko & Kallman 1994). This method consists 
in keeping the gas pressure constant during the iterations 
for convergence of the energy balance equation; this proce- 
dure is correct, but leads to a situation of thermal instabil- 
ity, as described in the previous sections. Photoionization 
codes using this computational method (e.g. XSTAR: 
Kallman & Krolik 1995; Kallman & Bautista 2001) must 
then choose arbitrarily one of the stable solutions in the 
multi-branch region. The second possibility is to keep the 
density constant during the iterations for convergence of 
the energy balance equation, allowing to obtain the tem- 
perature value in the layer; this computational method 
provides a unique solution, even if it is an approximate 
one, intermediate between the "cold" and the "hot" 
solutions present in the multi-branch region. This ap- 
proach has been used in studies of the vertical structure 
of accretion discs (Raymond 1993; Shimura, Mineshige 
& Takahara 1995; Madej & Rozahska 2000; Kawaguchi, 
Shimura & Mineshige 2001; Ballantyne et al. 2001). 

Such a computational method is implemented in our 
photoionization code TITAN, described in the following 
section. It has been used during the past years, with good 
results; models computed with the "intermediate" tem- 
perature solution have been applied, for instance, to stud- 
ies of the vertical disc structure with hydrostatic equilib- 
rium (Rozahska et al. 2002), and to total pressure equi- 
librium media (Rozahska et al. 2002; Dumont et al. 2002; 
Rozahska et al. 2006; Gongalves et al. 2006; Chevallier et 
al. 2006a). Since then, we have implemented a new algo- 
rithm in the TITAN code; here, we report on this new 
addition and its applications. In summary, the new algo- 
rithm allows to choose between the "hot" or the "cold" 
stable solutions in the multi-branch region, and thereof 
to compute a fully consistent model using the chosen so- 
lution. For the first time, it is thus possible to compare 
models issued from the true stable solutions with the ap- 



4 This "intermediate" solution should not be identified with 
the unstable solution (for which the derivative of L is negative); 
it is simply a numerical alternative, intermediate between the 
stable solutions. 



proximate solution model and to estimate the degree of 
uncertainty of such an "intermediate" solution. Our im- 
proved version of the TITAN code will allow to better 
understand the behavior of both kind of computational 
schemes and to estimate their influence on the description 
of the ionized gas structure and on the modelled spectra 
emitted and absorbed by the medium. 

In Sect. [2] we give a brief description of the TITAN 
code and discuss the different computational methods 
used in more detail; the models are described in Sect. [3l 
Section[4]contains a comparative study between the "hot" , 
"cold", and "intermediate" solutions. In Sect. [5] we discuss 
some observational implications, and in Sect. [6] we sum- 
marize our conclusions. 

2. Computational issues 

2.1. The TITAN code 

All our models, summarized in Table 1, were computed 
using the TITAN code. TITAN is a photoionization- 
transfer code developed by our team to correctly model 
optically thick (Thomson optical depth up to several tens) 
ionized media; it can equally be applied to thinner me- 
dia (Thomson depth ~ 0.001 to 0.1, e.g., Collin et al. 
2004; Goncalves et al. 2006). The code includes all rel- 
evant physical processes (e.g., photoionization, radiative 
and dielectronic recombination, fluorescence and Auger 
processes, collisional ionization, radiative and collisional 
excitation/de-excitation, etc.) and all induced processes. 
It solves the ionization equilibrium of all the ion species 
of each element, the thermal equilibrium, the statistical 
equilibrium of all the levels of each ion, and the transfer of 
both the lines and the continuum. It gives as output the 
ionization, density, and temperature structures, as well as 
the reflected, outward, and transmitted spectra. The en- 
ergy balance is ensured locally with a precision of 0.01%, 
and globally with a precision of 1%. 

TITAN is based on an idea initially depicted in Collin- 
Souffrin & Dumont (1986); its conception took several 
years, the code being described for the first time in 
Dumont, Abrassard & Collin (2000). It has been perma- 
nently upgraded since then; some of the improvements 
are described in, e.g., Dumont & Collin (2001), Dumont 
et al. (2002), and Chevallier et al. (2006a). In particular, 
our code uses the Accelerated Lambda Iteration (ALI) 
method, which allows for the exact treatement of the 
transfer of both the continuum and the lines (see Dumont 
et al. 2003 for a description of the ALI method in the 
modelling of the X-ray spectra of AGN and X-ray bina- 
ries). This is a major improvement over other photoioniza- 
tion codes such as Cloudy (Ferland et al. 1998), XSTAR 
(Kallman & Krolik 1995; Kallman & Bautista 2001), or 
ION (Netzer 1993, 1996) which use, at least for the lines, 
an integral formalism called the "escape probability ap- 
proximation" . While ALI computes very precisely line and 

5 Our atomic data include ~ 10 lines from ions and atoms 
of H, He, C, N, O, Ne, Mg, Si, S, and Fe. 
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continuum fluxes, treating them in a consistent way, in 
approximate methods the computation of the absorption 
and emission lines is uncoupled. Furthermore, ALI enables 
a multi-direction utilisation of the code, allowing for any 
illumination angle (including the normal direction) and 
any direction of the outward and reflected emission. This 
fully operational version of the code offers the possibility 
to model an ionized medium in constant density, constant 
gas pressure, or constant total (i.e. gas plus radiation) 
pressure. It has recently been used to study the soft X-ray 
spectra of AGN (Chevallier et al. 2006a) and applied to 
model real data as, for instance, the Warm Absorber in 
NGC 3783 with a single medium in total pressure equi- 
librium (Gongalves et al. 2006), or the spectra of bright 
ultraluminous X-ray sources (Gongalves & Soria 2006). 

2.2. Different computational methods used 

We recall that the set of equations describing the phys- 
ical state of the gas (i.e. the local balance between the 
ionization and recombination processes, excitations and 
de-excitations, as well as the local energy balance) is com- 
puted at each depth in the gas slab. In TITAN's previous 
computational scheme, the thermal energy balance equa- 
tion was solved for a given density at a given depth (that 
of the previous iteration); this method provides a unique 
solution for the temperature. The model corresponding to 
this scheme will be called "intermediate" model hereafter. 



Model 


£ = L/nuR' 


Nn 


Thermal instability 


name 


(erg cm s" 1 ) 


(cm" 2 ) 


solution 


HI1_C 


10000 


210 23 


cold 


HI1_H 


10000 


210 23 


hot 


HI1J 


10000 


210 23 


intermediate 


MI1_C 


1000 


210 zd 


cold 


MI1J3 


1000 


210 23 


hot 


MI1J 


1000 


210 23 


intermediate 


MI2_C 


1000 


2.5 10 23 


cold 


MI2JH 


1000 


2.5 10 23 


hot 


MI 2 J 


1000 


2.5 10 23 


intermediate 


MI 3_C 


1000 


310 23 


cold 


MI3_H 


1000 


310 23 


hot 


MI 3 J 


1000 


310 23 


intermediate 


LI1_C 


300 


510 22 


cold 


LI1_H 


300 


510 22 


hot 


LIU 


300 


510 22 


intermediate 


LI2_C 


300 


8.5 10 22 


cold 


LI2_H 


300 


8.5 10 22 


hot 


LI2_I 


300 


8.5 10 22 


intermediate 



Table 1. Models computed with TITAN in the frame- 
work of this study. In the first column, the model 
name regroups information on the ionization (H: high, 
M: medium, L: low) and on the chosen computational so- 
lution (C: cold, H: hot, I: intermediate). The ionization 
parameter and column density values are given in Cols. 
2 and 3, respectively. Col. 4 lists the thermal instability 
solution chosen. 
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10 100 1000 10000 

hv (eV) 

Fig. 1. This figure displays the spectral distribution of 
the mean continuum intensity, computed at the illumi- 
nated surface of the slab, for three "intermediate" models: 
MI 1_I (dotted line), MI 2 J (solid line), and MI 3_I (dashed 
line). For the sake of clarity, the spectral lines have been 
supressed from this figure, where only the continuum is 
shown; this displays important discontinuities, which in- 
crease with the slab thickness. 



In the new version of TITAN, the energy balance 
equation is solved for a given total pressure, assumed to 
be constant along the slab. In some regions of the cloud, 
this computational scheme leads to two stable solutions, 
described by two separate models. In one of the models, 
the "hot" solution is systematically chosen for all layers, 
while in the other model the "cold" solution is system- 
atically preferred. These models will be called hereafter 
"hot" model and "cold" model, respectively. They actu- 
ally represent two extreme results corresponding to a given 
gas composition and photoionizing flux, being compatible 
with the two stable solutions. It is important to stress that 
these models differ not only in the layers where multiple 
solutions are possible, but all along the gas slab, this be- 
cause the entire radiation field will be modified in a thick 
medium. 

This, and other behaviours of a slab of optically thick 
gas in constant pressure equilibrium, will be discussed in 
the following sections. We have chosen to work with a few 
models encompassing the range of conditions valid for the 
Warm Absorber in AGN; they are given in Table 1. The 
same study could, of course, be performed for a different 
set of values; this would not change our main results and 
conclusions. 

3. The models 

All the models presented here were treated in a 1-D plane- 
parallel geometry, as slabs of gas illuminated from one side 
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Fig. 2. This figure displays some results obtained for an "intermediate" model with £ = 1000 and A^h = 2 10 cm 
(Mil J). On the left-hand panel we give the spectral distribution of the mean continuum for two layers located at 
different depths in the gas slab; these layers correspond to a column density value of 1.66 10 23 (dashed line) and 
1.90 10 23 (solid line). On the right-hand panel, we display the S-curves giving T versus £/T for the same two layers 
(the same line-codes apply); the dots represent the equilibrium temperature found for each layer. 



by a radiation field concentrated in a very small, pencil- 
like, shape centered on the normal direction. The resul- 
tant spectrum, reprocessed by the gas, is a combination 
of "reflection" from the illuminated side of the medium 
(not a real reflection, as it includes atomic and Compton 
reprocessing), "outward emission" (coming from the non- 
illuminated side of the medium), and a "transmitted" frac- 
tion of the incident ionizing continuum. The absorption 
spectrum (which corresponds to the outward emission in 
the normal direction) depends strongly on the column den- 
sity, since the cool absorbing layers are located at the back 
side of the gas slab. The relative contribution of each com- 
ponent to the observed spectrum depends on parameters 
such as the size, the density, and the geometry of the ion- 
ized medium. 

Our models are globally parameterized by the column 
density A^h (in units of cm -2 ), the hydrogen number den- 
sity at the illuminated side of the slab «h (in units of 
cm -3 ), the incident flux and SED, and the ionization pa- 
rameter; this parameter is defined differently among au- 
thors. In this article, we use the £ form of the ionization 
parameter, where 

£ = - (in units of erg cm s -1 ) (3) 

nuR 2 

in which R is the distance between the radiating source 
and the photoionized medium, and L is the source's bolo- 
metric luminosity (in erg s -1 ); in the following, we inte- 
grate L over the range 10— 10 5 eV, but some authors prefer 
to give £ with L integrated over the 1—1000 Ryd region, as 



used by XSTAR. Appropriate conversions should thus be 
applied if using different forms of the ionization parame- 
ter, and before any numerical comparison. Note that £, as 
well as nfj, are quantities defined at the illuminated sur- 
face of the slab. We recall that in constant pressure mod- 
els the hydrogen number density varies along the slab, 
displaying a profile varying inversely with the tempera- 
ture; also the spectral distribution of the ionizing radiation 
changes across the slab for media in pressure equilibrium, 
as those discussed in this paper. 

All our models were computed under total (i.e., gas 
plus radiation) pressure equilibrium; we have neglected 
all other forms of pressure (e.g., turbulent, magnetic, ...) 
contributing to the total pressure. The gas slabs were as- 
sumed to have cosmic abundances (Allen 1973); the hy- 
drogen numerical density at the illuminated surface was 
set to 10 7 cm -3 . Note that n# is a minor parameter in 
this study, the overall spectrum being only proportional 
for values varying from 10 7 to 10 12 cm -3 (this is usually 
the case, provided the ionizing spectrum is relativelly flat) ; 
only the relative intensities of the forbidden and permit- 
ted emission lines (e.g., from He- like ions like O vn) are 
sensitive to changes in the density. 

We have simplified our models by assuming a null tur- 
bulent velocity component; micro-turbulence has a small 
influence on the thermal and ionization structure for 
large column densities, where the cooling is dominated by 
bound-free transitions. The thermal and ionization struc- 
tures are mainly determined by the ionization parameter 
and the shape of the ionizing continuum. All the models in 
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Fig. 3. On the left-hand panel we give the temperature profile along the slab (i.e., versus the column density), for 
the model MI 2 (£ = 1000, Nn = 2.5 10 23 ) computed with the "cold" solution. The dots represent gas at different 
depths in the slab, for which we have plotted the corresponding cooling curves on the right-hand panel. We observe 
that the net cooling function C has one single solution for the layers labelled 0, 1 and 4 in the slab, while there are 
three possible solutions for the layers labelled 2 and 3; from these, only the "cold" solution was selected to compute 
the gas structure and spectra. 



this paper assume an incident continuum described by a 
power-law of photon index T — 2, covering the 10— 10 5 eV 
energy range. 

4. "Hot", "cold", and "intermediate" solutions 

In this section, we illustrate the differences between the 
previously discussed computational methods by compar- 
ing the results obtained with the "hot" , "cold" , and "in- 
termediate" models. 

4.1. Behaviour at the surface of the gas slab 

It is often assumed that the thermal instability problem 
in thick media is exactly the same as in thin media. This 
is not true, for mainly two reasons. First, the spectral 
distribution of the mean intensity J„ at the illuminated 
surface of the slab is different from the incident spectrum, 
as it equally contains a "returning" radiation component 
emitted by the slab itself. As an illustration, Fig. [1] shows 
the spectral distribution of the mean continuum intensity 
at the slab surface for three "intermediate" models with 
£ = 1000 and 3 different values of the column density. For 
the sake of clarity, the spectral lines were suppressed from 



the figure, and only the continuum is shown. We observe 
that the spectrum at the surface contains strong discon- 
tinuities, whose amplitude increases with the slab thick- 
ness. More generally, the intensity of the whole spectrum 
increases with the slab thickness (an expected behaviour 
since the radiation emitted by a thicker slab should be 
larger). Second, the spectral distribution changes as the 
radiation progresses inside the medium. As a consequence, 
the shape of the S-curve also changes, and instead of trav- 
eling along a given S-curve as P r ad/-Pgas decreases, the 
temperature follows successively different curves. This be- 
haviour is illustrated by the MI U model in Fig. [2j which 
shows the spectral distribution of the mean continuum for 
two layers located at different depths in the gas slab (left- 
hand panel) and the corresponding curves giving T versus 
£/T (right-hand panel), which is the same as represent- 
ing the usual S-shape curve. We can see that the S-curves 
are different for the two represented layers and, in par- 
ticular, that the curve corresponding to the deeper layer 
(solid line) displays a larger multi-branch region; this is 
due to its larger absorption trough. The dots represent 
the equilibrium temperature found for each layer. 
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Fig. 4. Same as in Fig. 3, except the model has been computed with the "hot" solution. In the right-hand panel, layers 
labelled 2 and 3 display three possible solutions, from which the "hot" one has been systematically choosen; the layers 
labelled 0, 1, and 4 display a single solution. 




Fig. 5. Same as in Figs. 3 and 4, except the model has been computed with the "intermediate" solution. In this case, 
the cooling curves in the right-hand panel always provide a unique solution, independently of the depth of the layer 
considered (dots on the left-hand panel). 



4.2. Comparison between cooling curves 

Figures [31 [U and [5] concern a slab of gas with total col- 
umn density equal to 2.5 1 23 cm -2 , ionized by an inci- 
dent continuum with £ = 1000 (MI 2 model). These fig- 



ures display, for respectively the "cold" , "hot" , and "inter- 
mediate" solution models, the temperature profile versus 
the column density (left-hand panel) and the net cooling 
function (C = A — T) versus the temperature, at different 
depths in the gas slab (right-hand panel); the various lay- 



8 



A. C. Gongalves et al.: Thermal instability in X-ray photoionized media in AGN 




le+22 2e+22 3e+22 4e+22 5e+22 6e+22 7e+22 8e+22 
N„(cm 2 ) 




I.Se+23 2e+23 2.5e+23 



le+23 1.5e+23 
N,,(cm 2 > 



Fig. 6. Temperature profiles versus the column density 
for a set of models with different £ and Nh. On the top 
panel, the LI 1 models (dashed line: "hot" solution; dotted 
line: "cold" solution; solid line: "intermediate" solution); 
for comparison, we also plotted the LI 2_L model (dashed- 
dotted line). On the middle panel, several MI models for 
different values of the column density (dashed lines: "hot" 
solution models; dotted lines: "cold" solution models; solid 
lines: "intermediate" solution models). Finally, on the bot- 
tom panel, the HI1 models (the same line-codes apply). 



ers represented here have been marked and labelled in the 
left-hand panel. By comparing both plots we can verify 
that C cancels for the equilibrium temperatures observed 
in the gas structure. 

We recall that the cooling curves of the "hot" and 
"cold" models are determined for a given total pressure, 



while those of the "intermediate" models are computed 
for a given density. We can check that the isobaric cool- 
ing curves of the "cold" and "hot" models (Figs. [3] and 
IH respectively) display three solutions for C = 0, while 
the isodensity cooling curves of the "intermediate" model 
(Fig. [5]) always show a unique solution. We stress that 
it is impossible to compare the cooling curves in exactly 
the same conditions, as the models computed with these 
three solutions differ considerably. Indeed, the radiation 
field is completely different in the three models, even at 
the illuminated surface layer; therefore, the isobaric cool- 
ing curves differ strongly, even for a layer located at a 
similar depth in the gas slab. 

Let us examine first the cooling curves corresponding 
to the "cold" model (Fig. [3]). We can see that, at the 
illuminated surface and at position 1, there is a unique 
solution, while at positions 2 and 3 there are three possi- 
ble solutions, from which the coldest one is chosen; near 
the back surface of the slab, at position 4, the solution is 
again unique. Similarly, for the "hot" model (Fig. 0J, the 
solution is unique at the surface and at position 1, but 
at positions 2 and 3 the curve has three solutions, from 
which the hotter one is chosen; then, at position 4, again 
the solution is unique. However, the multi-branch region 
is different in the case of the "hot" or the "cold" model; 
this shows how much these two models differ. Finally, in 
the case of the "intermediate" model (Fig. [5]), all solu- 
tions are unique. Given these differences, we would expect 
that the previous computational scheme used by TITAN, 
with its isodensity cooling curves, would be completely 
wrong; this is not the case, as we will see later on. Indeed, 
such a method provides temperature structures and spec- 
tra which are intermediate between those of the "cold" 
and the "hot" models. 

4.3. Temperature profiles 

Figure [6] displays the temperature profiles for the "cold" , 
"hot" , and "intermediate" solutions, for a set of LI (£ = 
300), MI (£ = 1000), and HI (£ = 10000) models. The "in- 
termediate" solutions are represented by solid lines (on the 
top panel, the LI 2_I model was represented by a dashed- 
dotted line for the sake of clarity) , the "cold" solutions by 
dotted lines, and the "hot" solutions by dashed lines. By 
analysing this figure, two results appear clearly: 

— First, the differences between the "hot", "cold", and 
"intermediate" models increase as £ decreases; 

— Second, the temperature always decreases strongly 
near the non-illuminated surface of the slab, whatever 
the column density. 

In the following, we will present some clues allowing 
to understand these behaviors. 



4.3.1. Influence of the ionization parameter 

For the higher ionization models (£ = 10000), the "hot", 
"cold" , and "intermediate" solutions provide very similar 
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Fig. 7. Temperature profile versus the column density for a model with £ = 300 and N-r = 5 10 22 (left-hand panel). 
We have marked several layers, corresponding to different positions in the gas slab, labelled from (at the illuminated 
surface) to 4 (deeper in the slab). On the right-hand panel, we give the corresponding cooling curves, which display 
one (position 4) or more solutions. It is interesting to note that, at such low ionizations, the surface layer is already in 
the multi-solutions regime, displaying 3 possible solutions from which the coldest one has been chosen. 
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Fig. 8. This figure shows various pressure quantities computed with the "cold" (thick lines), "hot" (thin lines), and 
"intermediate" (thinner lines) solutions for the models LI 1 (left-hand panel), MI 1 (middle panel), and HI 1 (right-hand 
panel). For the sake of clarity, the "intermediate" solution values are only shown for one of the models (MI 1, middle 
panel). In all the panels, P ra d is given in solid lines, P gas in dashed-dotted lines, P ra d (lines) in dotted lines, and the 
P rac j/P gas ratio in dashed lines; also, the left-hand y-axis (10 -4 — 1) corresponds to the P ra< a, P gaS j & n d P rac i (lines) 
values, while the right-hand y-axis (10 -2 — 10 2 ) corresponds to the P ra d/-Pgas ratio. 



results (cf. Fig. [HI bottom panel). On the contrary, for the 
lower ionization models (£ = 300) the results are quite 
different, the first layer at the illuminated side of the slab 
being already in the multi-solutions regime (cf. Figure El 
top panel). In this case, since the "hot" and "cold" so- 
lutions correspond to very different equilibrium tempera- 
tures, there is also a large difference on the whole medium 



structure. This is illustrated by Fig. \7\ which is similar to 
figures 3, 4 and 5, except we now address the case of the 
LI 1 "cold" solution model. The figure shows that, already 
at the illuminated surface of the slab (labelled on the 
left-hand panel), the net cooling curve cancels for three 
values of the temperature; this multi-solutions regime lasts 
until the deeper position labelled 4 on the left-hand panel. 
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Fig. 9. Figure showing the variation of P ra d, P- a d (lines), 
and P gas versus the column density, for 3 "intermediate" 
models of different thickness (MI 1, MI 2, and MI 3). P rad 
is represented by solid lines, P gas by dashed lines, and 
P ra d (lines) by dashed-dotted lines. The variation of P ra d 
along the gas slab is very similar for the three different 
thicknesses, until the first temperature jump occurs. We 
also observe that near the back side of the gas slab, P ra d 
can become very important, owing to the contribution 
from the lines. The thin dotted vertical lines represent lay- 
ers at different depths in the slab, corresponding namely to 
1 10 23 cm -2 , 1.5 I0 23 cm~ 2 , and to the first temperature 
jump in each model (1.82 10 23 , 2.16 10 23 , and 2.43 10 23 
cm~ 2 for the Mil, MI 2, and MI 3 model, respectively). 



It is possible to address this situation and try to find an 
adequate explanation for the behaviour of the LI 1 model 
by looking into what's happening with the 5 ionization 
parameter (or, which is similar, studying the P ra d/Pga,s 
ratio). Since the multi-solutions regime occurs for a given 
range in S, one expects that the associated region in the 
slab would correspond to layers where P ra d/P g as is of the 
order of a few units. This is indeed the shown 
in Fig. [U This figure displays the variation of P ra d, Pgas, 
and P ra d/P ga s along the slab, for the "hot" and "cold" 
solutions of a set of LI, MI, and HI models. The radia- 
tion pressure due to the lines, P ra d (lines), is also shown, 
as this component can be quite important near the back 
side of the slab; note that P ra d, also represented in the fig- 
ure, already includes the contribution from P rad (lines). In 
this figure, the "cold" , "hot" , and "intermediate" models 
are represented by thick (green, black, violet), thin (blue, 
orange, magenta), and very thin lines (cyan, yellow, red), 
respectively. P rad is given in solid lines, P gas in dashed- 
dotted lines, P rad (lines) in dotted lines, and the P ra d/Pg as 
ratio in dashed lines. The left-hand y-axis (I0~ 4 — I) cor- 
responds to the Pad, Pgas, and P ra d (lines) values, while 




1.5e+23 2e+23 2.5e+23 3e+23 



Fig. 10. Figure showing the P-ad/Pgas ratio versus the 
column density, for the 3 "intermediate" models discussed 
in the previous figure. This ratio decreases as we penetrate 
in the gas slab until it reaches a value ~2.5, coinciding 
with a temperature jump. The figure shows that this ratio 
decreases more rapidly for the thinner, than for the thicker 
gas slab. The thin dotted vertical lines represent layers at 
different depths in the slab (same as in Fig. 9). 



the right-hand y-axis (I0~ 2 — I0 2 ) corresponds to the 
Prad/Pgas ratio. For the sake of clarity, the "intermediate" 
solution values (thinner lines: cyan, yellow, red) are only 
shown for one of the models (MI I, in the middle panel). 

Observing Fig.[8l we first notice that, in all cases, the 
pressure is entirely dominated by P ra d, and that P ra d/Pgas 
at the illuminated surface increases with £, as expected. 
The P-ad/Pgas ratio decreases accross the slab from the 
illuminated surface to the back side as the radiation is ab- 
sorbed, the slab entering the multi-solutions regime when 
P-ad/Pgas is of the order of 2.5; this is true whatever the 
pressure ratio at the illuminated side. It is easy to see that 
such a value was promply attained in the case of the LI I 
model discussed before; as a consequence, we were able 
to observe a "cold" and a "hot" solution already at the 
surface of the slab. 

It is also interesting to note that, in spite of a big dif- 
ference in temperature, the gas pressure at the surface is 
almost the same for the "cold" and "hot" models. This 
behaviour is actually due to an automatic adjustment of 
the density, which is equal to 5 I0 7 cm~ 3 in the first layer 
of the "cold" model, instead of the smaller, initial value of 
I0 7 cm~ 3 . The figure also shows that, for all values of the 
ionization parameter, the pressures in the three models 
("cold", "hot", and "intermediate") are quite similar, ex- 
cept for P a d (lines) near the back side of the slab. Indeed, 
the contribution from the lines is always less important 
in a hot medium than in a cold medium, which can emit 
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Fig. 11. Figure showing the medium intensity spectral distribution for two MI "intermediate" models with different 
column densities: on the left-hand panel, the thiner (shorter) model with A^h = 2 10 23 ; on the right-hand panel, the 
thicker (longer) model with A^h = 3 1 23 . Both panels display the spectral distribution corresponding to 4 layers, from 
the illuminated surface of the slab (top curves) to the layer closer to the temperature jump (bottom curves). 



intense ultraviolet lines. Therefore, the smaller the ion- 
ization parameter, the more important are the lines. It is 
thus easy to understand how P ra d (lines) can even domi- 
nate the radiation pressure near the back side of the the 
LI1_C model (cf. Fig. HI left-hand panel). 

Note finally that the sharp temperature drops, ob- 
served in the left-hand panel of figures 3, 4, 5 and 7, do 
not occur exactly at the same value of P ra d/Pgas for all the 
models; this is due to the shape of the S-curve, and to the 
different spectral distributions inside the slab, which de- 
pend not only on the type of model ("cold" or "hot"), but 
also on the ionization parameter and on the total tickness 
of the medium. 

4.3.2. Influence of the medium thickness 

As already noted by Rozariska et al. (2006) and Chevallier 
et al. (2006a), the thickness of a pressure equilibrium 
medium cannot exceed a maximum value for a given ion- 
ization parameter. This is explained by the fact that, when 
the temperature reaches low values (of the order of 10 4 K), 
the radiation pressure becomes dominated by the spectral 
lines, and therefore the gas enters in a new multi-solutions 
regime. When the temperature drops to the cold values 
corresponding to a molecular gas, we reach a limit and 
are unable to provide a solution, as our code cannot han- 
dle molecular gas. 

If the maximum thickness of a medium in pressure 
equilibrium is a known issue, there are two other puzzling 
problems related to the thickness of the slab. The first is 
related to the fact that the temperature always decreases 



rapidly near the non-illuminated side of the gas slab. The 
second concerns the reason why the slab always enters 
a multi-solutions regime just before the slab ends, even 
when the imposed column density is much smaller than 
the possible maximum value. 

We are going to address these issues using the "inter- 
mediate" solution models as an illustration. Let us con- 
sider a set of models with £ = 1000. Figure M displays 
the variation of P ra d, P-ad (lines), and P gas along the slab 
for three values of the column density. P ra d is represented 
by solid lines, P gas by dashed lines, and P rac j (lines) by 
dashed-dotted lines. We can see that the dependence of 
P-ad — and therefore also of P gas , since the sum of P ra( j 
and Pg as is constant — on the gas position inside the slab 
is almost the same for the three different thicknesses until 
the first temperature jump occurs. For the sake of sim- 
plicity, we do not show here P-ad/Pgas, but it is clear that 
this ratio displays always approximately the same value 
(~ 2.5) at the temperature jump (the position of the tem- 
perature jump can be seen in Fig. [6l middle panel). Such 
temperature changes are accompanied by a strong increase 
in P ra d (lines). The thicker the slab, the larger is the cold 
zone, and the stronger the increase in P rat j (lines). One 
can even observe an inversion of the gas pressure varia- 
tion near the back surface of the slab, for the thickest case 
(MI 3 model). Very close to the back surface, the layers 
are optically thin for the continuum and the lines, which 
can thus escape from the medium; this induces a rapid de- 
crease in the radiation pressure, as observed in the figure; 
it also leads to a strong cooling, causing a rapid decrease 



12 



A. C. Gongalves et al.: Thermal instability in X-ray photoionized media in AGN 



le+08 



le+07 




le+08 



100 1000 
hv (eV) 



10000 



> 



le+07 




1200 



hv (eV) 



le+08 i 



le+07 



> le+06 



le+05 



i 1 1 in 1 — i'ii 'in 1 1 1 — i — i i i 1 1 




_i i 



_i 



_i 



10 



100 



1000 



10000 



hv (eV) 



le+08 



le+07 



> le+06 



le+0: 




1200 



hv (eV) 



Fig. 12. This figure gives the pure absorption (top panels) and outward emission spectra (bottom panels) for the 
MI 1 model. The "cold" solution model is represented by the lower blue lines, the "hot" model by upper magenta 
lines and the "intermediate" model by the middle black lines; the spectra corresponding to the half-sum of the "hot" 
and "cold" models are also plotted for comparison (green dotted lines). The orange dashed line corresponds to the 
incident ionizing continuum. The spectra are represented with a resolution of 300; the right-hand panel foccuses on 
the 700—1200 eV energy range for details, while the left-hand panel shows a larger band covering the 10— 10 4 eV 
range. Notice the excellent agreement between the "intermediate" and half-sum spectra over the entire energy range. 



of the temperature near the back surface of the slab. These 
coupled phenomena give an answer to the first problem. 

The second issue is more subtile. From Fig. [9] we know 
that P ra a decreases inside the slab (except near the back 
side, owing to the lines, as explained). Since P ra d is slightly 
larger in thicker slabs, the traversed length before reach- 
ing the ratio P ra d/ Pgas required for the temperature jump, 
is larger for thicker slabs. However, this is not sufficient to 



explain why the temperature transition is located closer 
to the illuminated surface in thiner slabs, or the difference 
between the transition positions would not be so large in 
media with different thicknesses. Additional explanation 
can be found in Fig. [TUl which compares the P r ad/-fgas 
ratio versus the position in the slab (given in terms of col- 
umn density) for the same three slab thicknesses discussed 
in the previous figure; we can see that this ratio decreases 
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Fig. 13. Same as in the previous figure, only for the lower ionization model LI 1. In this case, the agreement between 
the "intermediate" and half-sum spectra is not as good as for the MI 1 model. The differences arc less important for 
the outward emission in the X-ray range. 



more rapidly for the thinner (i.e. shorter), than for the 
thicker (longer) gas slab. 

The explanation for this behaviour can be found in the 
comparison between the spectral distribution of the first 
layer, and of a layer located close to the place where the 
jump in temperature occurs. Figure QT] displays the spec- 
tral distribution of the mean intensity (again, the lines 
were suppessed for clarity), for two MI, "intermediate" 
models with total column densities 2 10 23 cm -2 (left-hand 
panel) and 3 10 23 cm -2 (righ-hand panel), at different po- 
sitions in the slab (from top to bottom: at the illuminated 
surface, at 1 10 23 cm -2 , 1.5 10 23 cm -2 , and at the first 



temperature jump — 1.82 10 23 for MI 1 and 2.43 1 23 for 
MI 3); these layers are identified in Fig.QJJlby thin vertical 
dotted lines. We can observe that the spectral distribution 
of the continuum varies differently in the two media: at the 
same location in the slab, the spectrum is less absorbed 
in the thicker slab (righ-hand panel), owing to the more 
important re-emission component. In other words, the 
spectral distribution of the thiner, shorter slab (left-hand 
panel) has a stronger trough, and reaches more rapidly the 
multi-solutions regime, as explained above (cf. Fig. [J). 
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5. Observational implications 

5.1. Outward emitted and absorbed spectra 

Unfortunately, it is impossible to know what solution 
the plasma will adopt when attaining the multi-solutions 
regime. For instance, it can oscillate between the hot and 
cold solutions, it can fragment into hot and cold clumps 
which will coexist together, or it can take the form of 
a hot, dilute medium confining cold, denser clumps (the 
hot and cold media sharing the same pressure). We recall 
that the "cold" and "hot" models computed by TITAN 
can be interpreted as "extreme" models, as they adopt re- 
spectively the "cold" or the "hot" solutions all along the 
slab. The spectra emitted or absorbed by a given ionized 
medium, consisting of a mixture of gas in the hot and cold 
phases, should thus be intermediate between those result- 
ing from the pure "cold" and "hot" models. The differ- 
ences between the emitted and absorbed spectra obtained 
with the stable solutions can thus provide an indication 
for the maximum "error bars" associated to the spectra 
computed with the "intermediate" solution. 

As an example, Figs. [T2l and [T3l give the pure absorp- 
tion (top panels) and outward emission (bottom panels) 
spectra resulting from the "cold" , "hot" , and "interme- 
diate" (black, middle curves) solutions for the MI 2 and 
LI 1 models, respectively (see Fig. [5] for the corresponding 
temperature profiles). The spectra are represented with a 
resolution of 300 and have been displayed in two energy 
ranges: 10— 10 4 eV (left-hand panels), with a zoom in the 
700—1200 eV region (right-hand panel) to more clearly 
observe the differences. The lower solid lines (blue) cor- 
respond to the "cold" models, the middle lines (black) to 
the "intermediate" models, and the upper lines (magenta) 
to the "hot" models. The incident ionizing continuum (or- 
ange dashed lines) and the spectrum corresponding to the 
half-sum of the "cold" and "hot" models (green dotted 
lines) are also plotted, for comparison. The outward emis- 
sion spectrum has been computed for an opening angle 
of 45 degrees which, according to the unified Scheme of 
AGN (Antonucci & Miller 1985), is in agreement with the 
expected geometry for the Warm Absorber in Seyfert 1 
nuclei, or the emissive region in Seyfert 2s. 

Investigating the MI 2 spectra in Fig. Q2] we can see 
that the "intermediate" model spectra (black solid lines, 
middle) is extraordinary close to the half-sum (green dot- 
ted lines) of the "cold" and "hot" spectra, though these 
ones are quite different. The differences amount at most 
to 1 or 2% for the outward emitted spectrum, and to 20% 
for the absorption spectrum. The agreement is not as good 
for the second case (LI 1 model, Fig. [13]), where the differ- 
ences between the "intermediate" model and the half-sum 
of the "cold" and "hot" models amount typically to a fac- 
tor 2 for the absorption spectrum. The agreement is much 
better for the outward emission in the X-ray range. This is 
due to the strong imprint of the cold layers near the back 
surface of the "cold" model, which absorb almost com- 
pletely the UV-X spectrum. On the contrary, the X-ray 



outward emission, being due mainly to the deep layers, is 
not very different for the "hot" and the "cold" models. 

We observe that the lines have similar intensities in the 
"cold" and "hot" spectra, and also in the "intermediate" 
and half-sum spectra; this happens in spite of important 
differences in the continuum intensity, in particular for the 
ionization edges. Such a behaviour is due to the fact that 
the main lines are saturated resonance transitions, while 
the ionization edges are optically thiner, and therefore al- 
most proportional to the thickness of the region containing 
the emitting ion (apart from the diffuse emission contri- 
bution). 

5.2. Ionization states 

Differences in the ionization state for "hot" and "cold" 
models are illustrated by Fig. [14] which displays the frac- 
tional ionization of two important elements — Si (top 
panels) and O (bottom panels) — for the "cold" (left- 
hand panels) and "hot" (right-hand panels) models al- 
ready depicted in Fig. [12] This figure clearly shows why 
the intensities of the ionization edges are different in the 
"cold" and "hot" models. We can see that the length of 
the O vii and O vni regions is about three times larger 
in the "cold" model than it is in the "hot" model, these 
differences being somewhat larger than the discrepancy in 
the corresponding spectra. 

5.3. Variability 

In the multi-solutions regime, the medium should be made 
of a mixture of gas in the hot and cold phases, and the 
relative proportion of those phases could be varying with 
time. The whole medium should thus oscillate between 
physical states located somewhere in the parameters' re- 
gion covered by the "hot" and "cold" models. The varia- 
tion timescale in such a medium is basically the thermal 
time in the multi-branch region. However, if the medium 
is subject to a change in temperature, this will not neces- 
sarily be immediately followed by a change in ionization 
inducing modifications in the main absorbing or emitting 
ions. It all depends on the relative values of the thermal 
time £th, the ionization time ti on and the recombination 
time t Icc . One thus expects two types of variability: a 
relativelly weak variability related only to a temperature 
variation, and a stronger variability related to a change of 
ionization state. If i; on or t TCC are smaller than tft, strong 
variations in the emitted/absorbed spectrum should be 
observed during the thermal time-scale, with the main 
emission/absorption lines and ionization edges being re- 
placed by lines and edges of other ions. If ti on and t rec 
are larger than tth, the spectrum should display smaller 
changes, like a simple variation in the relative intensity of 
the spectral features. 

Let us compare these three timescales. The ther- 
mal time t t h is of the order of kT/nA. In the multi- 
solutions regime, the cooling and the heating are domi- 
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nated by atomic processes, and A is of the order of 10 23 
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where ri\2 is the density expressed in 10 12 cm~ 3 , A23 
is the cooling function in 10 -23 erg cm 3 s _1 , and T5 is the 
temperature in 10 5 K. 

For a given ion, t rcc is equal to t ICC = l/(n c ai on ) where 
n e is the electron numerical density (roughly equal to the 
total density n), and ai on is the recombination coefficient 
of the given ion, including dielectronic recombinations. 
This coefficient is almost independent from the radiation 
flux, with a typical value of 10 -11 cm 3 s _1 at the tempera- 
ture of the multi-branch region for a heavy, highly ionized 
element, so t rec writes: 



0.1 



ni2 an 



(in units of s), 



(5) 



where the recombination coefficient an is expressed 
in 10~ n cm 3 s -1 . Thus, the recombination time t TCC of a 
heavy element in a highly ionized state is smaller than, or 
of the order of, the thermal time i t h- One deduces that a 
thermal variation induces almost immediately a change in 
the ionization state. 

Ionization equilibrium is reached after a time equal 
to the maximum time-scales given by U on and t rcc . The 
ionization time U on is defined as: 



i-i/(max) j ' 

/ 4-7T — — a v dv 

L(ion) hU 



(G) 



where <7„ is the frequency-dependent photoionization 
cross-section of the given ion. As we have seen, J v varies 
across the slab, but we can make a very rough estima- 
tion of ij on by assuming that J„ is of the order of the 
incident flux, F v /2n. This could lead to a large overesti- 
mation of the ionization rate, and therefore to an under- 
estimation of the ionization time, as J„ is actually smaller 
than F v /2-k in the deep layers, at the back of the slab (cf. 
the outward emitted spectrum in Figs. [T2l and [T3]) . With 
our power law incident spectrum described by F v oc v^ 1 
in the 10— 10 -5 eV range, and assuming that the cross- 
section is proportional to ^~ 3 , one findto 



3 IP' 5 , .. , , 
(in units of s), 

£3 ni2 CT19 



(7) 



where £3 is the ionization parameter £ expressed in 
10 3 erg cm s _1 and CT19 is the cross-section at the ionization 
edge of the given ion expressed in 10 -19 cm~ 2 . 

In conclusion, even if it may be underestimated by two 
orders of magnitude, t[ on is smaller than t lec , and both 
time-scales are smaller than i t h! therefore, if the temper- 
ature changes locally in the multi-branch region, a new 
ionization equilibrium settles in immediately. 



To derive this equation, we have made use of the definition 
of £ given in Eq. [3] 



The problem is more complex in practice. Indeed, in 
order to observe variations in the spectrum, the change of 
temperature should affect a large proportion of the multi- 
branch region. The mass fractions of the hot and cold 
phases are most probably of the same order, so this zone 
should be made of cold, dense clumps embedded in a hot- 
ter and more dilute medium. In this case, the evaporation 
time of the clumps should govern the structural changes. 
This time-scale is of the order of the cold phase dynamical 
time, if evaporation is not saturated. 

The dynamical time is roughly equal to AH/c s , where 
AH is the thickness of the multi-branch region, and c s is 
the sound velocity in the cold phase. The sound velocity 
includes the radiation pressure, but in a proportion which 
depends on the relation between the radiation field and 
the gas. Since the multi-branch region is not very opti- 
cally thick, one can neglect the radiation pressure, and 
simply assume c s = 2kT / inn, where toh is the proton 
mass. The thickness of the multi-branch region is a frac- 
tion /muiti of the total thickness of the slab, which itself 
depends on the density (since the column density is a given 
parameter). Expressing the total column density in 10 23 
cm -2 , one finds: 



idyn ~ 2 10' 



/multi Ni 



2:; 



(in units of s). 



(8) 



It is interesting to note that all timescales depend 
on the same way on the density, which is often an un- 
known parameter when dealing with observational data 
(e.g., Warm Absorber gas). If the timescales vary for dif- 
ferent locations in the slab (because the density value is 
not the same at different depths), their ratio does not 
change much. In particular, the ratio idynAth is very large 
(of the order of 10 3 if the thermal instability affects 10% of 
the slab), therefore the multi-branch region should contain 
local thermal perturbations out of pressure equilibrium. 

Chevallier et al. (2006b) made a detailed study on 
the viability of pressure equilibrium media in the pres- 
ence of variations of the incident radiation flux. The au- 
thors have shown that, though the global pressure equilib- 
rium is generally preserved (with the medium responding 
to the variations by rapidly adjusting its ionization and 
thermal equilibrium), supersonic velocities develop in the 
gas, owing to local over- or under-pressure states. The 
same phenomenon should occur in the multi-branch re- 
gion of our media, even in the absence of incident flux 
variations. To study this problem, a complete perturba- 
tional study should be performed around the equilibrium 
state of the gas slab; such a study is out of the scope of the 
present paper. For the time being, we can just conclude 
that large spectral fluctuations corresponding to the on- 
set of a "cold" or a "hot" solution could be observed in 
timescales of the order of the dynamical time, indepen- 
dently of any variation of the radiation field external to 
the medium. Moreover, a strong turbulence implying su- 
personic velocities leading to line broadening should per- 
manently exist in the multi-branch region of thick, strati- 
fied, pressure equilibrium media. 
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N H (cm" 2 ) N H (cm" 2 ) 

Fig. 14. Figure showing the fractional ionization of the Si ions (top panels) and O ions (bottom panels) for the 
same models as in Fig. Q21 The "cold" and "hot" models are displayed in the left-hand panels and right-hand panels, 
respectively. We observe that the regions producing the various ions have important size differences in the case of a 
"cold" , or of a "hot" model. 



5.4. Practical considerations 

An important point to take into account when choosing 
which computational method to apply, is that the full com- 
putation of the "hot" and "cold" models described here is 
extremely time-consuming; this is because the process is 
strongly unstable and requires thus more iterations than 
the isodensity scheme models. Given the results presented 
in this paper, it seems thus reasonable to use the simpler, 
isodensity scheme to compute constant pressure models 
or hydrostatic equilibrium models. This procedure is less 



ad-hoc than to choose arbitrarily between one of the pos- 
sible solutions resulting, in the end, in gas structures and 
emitted/absorbed spectra very close to what is expected 
from an ionized medium consisting in similar proportions 
of gas in the hot and cold phases. 



6. Conclusions 

1. We have addressed the thermal instability issue in 
the case of optically thick, stratified media, in total 
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pressure equilibrium. 

2. In order to do that, we have developped and im- 
plemented a new algorithm in the TITAN code; 
this algorithm, based on an isobaric computational 
method, allows to select the hot/cold stable solutions 
and thereof to compute a fully consistent photoionized 
model for each solution. 

3. We have chosen to work with a few models encom- 
passing the range of conditions valid for the Warm 
Absorber in AGN; our results can be applied to media 
in any pressure equilibrium conditions, e.g., constant 
gas pressure, constant total pressure, or hydrostatic 
pressure equilibrium. 

4. We have shown that a thick, stratified medium 
ionized by X-rays behaves differently from a thin 
ionized medium. This happens for mainly two reasons: 
first, the spectral distribution of the mean intensity 
at the illuminated surface of the slab is different 
from the incident spectrum, as it equally contains a 
"returning" radiation component emitted by the slab 
itself; second, the spectral distribution changes as 
the radiation progresses inside the medium and, as a 
consequence, the shape of the S-curve also changes. 
These effects depend on the thickess of the medium 
and on its ionization. 

5. This has observational implications in the emit- 
ted/absorbed spectra, ionization states, and variabil- 
ity. It is impossible to know what solution the plasma 
will adopt when attaining the multi-solutions regime: 
it can oscillate between the hot and cold solutions, 
it can fragment into hot and cold clumps which will 
coexist together, or it can take the form of a hot, dilute 
medium confining cold, denser clumps. Nevertheless, 
one expects the emitted/absorbed spectrum to be 
intermediate between those resulting from pure cold 
and hot models. 

6. We have compared the results obtained with models 
based on the pure hot /cold solutions, and models 
computed with an approximate, intermediate solution; 
we have demonstrated that the pure hot /cold models 
represent two extreme results corresponding to a 
given gas composition and photoionizing flux, being 
compatible with the two stable solutions. The three 
(hot, cold, and intermediate) models differ not only 
in the layers where multiple solutions are possible, 
but all along the gas slab, this because the entire 
radiation field suffers modifications while crossing a 
thick medium. 

7. The spectra emitted or absorbed by a given ionized 
medium, consisting of a mixture of gas in the hot 
and cold phases, should thus be intermediate between 
those resulting from the pure cold and hot models; 



17 

therefore, the intermediate model provides a good 
description of such a mixed-phase medium. The 
differences between the emitted and absorbed spectra 
obtained with the stable solutions provide an indica- 
tion for the maximum "error bars" associated to the 
spectra computed with the "intermediate" solution. 

8. The relative proportion of the hot and cold phases 
could vary with time. One expects two types of vari- 
ability: a relativelly weak variability related only to a 
temperature variation, and a stronger variability re- 
lated to a change of ionization state. Large spectral 
fluctuations corresponding to the onset of a cold/hot 
solution could be observed in timescales of the order of 
the dynamical time. Moreover, a strong turbulence im- 
plying supersonic velocities leading to line broadening 
should permanently exist in the multi-branch region of 
thick, stratified, pressure equilibrium media. 
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